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Abstract. We have extended our supercomputer-enabled Monte Carlo simulations 
of hopping transport in completely disordered 2D conductors to the case of substantial 
electron-electron Coulomb interaction. Such interaction may not only suppress the 
average value of hopping current, but also affect its fluctuations rather substantially. 
In particular, the spectral density Sj (/) of current fluctuations exhibits, at sufficiently 
low frequencies, a l//-like increase which approximately follows the Hooge scaling, 
even at vanishing temperature. At higher /, there is a crossover to a broad range of 
frequencies in which Si (/) is nearly constant, hence allowing characterization of the 
current noise by the effective Fano factor F = 5/ (/) /2e(/). For sufhciently large 
conductor samples and low temperatures, the Fano factor is suppressed below the 
Schottky value {F = 1), scaling with the length L of the conductor as F = {Lc/L)". 
The exponent a is significantly affected by the Coulomb interaction effects, changing 
from a ~ 0.76 ± 0.08 when such effects are negligible to virtually unity when they 
are substantial. The scaling parameter Lc, interpreted as the average percolation 
cluster length along the electric field direction, scales as oc £'-(0-98±o.08) .^^jjg^ 
Coulomb interaction effects are negligible and oc £'-(i-26±o.i5) -^^j^gjj such effects 
are substantial, in good agreement with estimates based on the theory of directed 
percolation. 



PACS numbers: 72.20.Ee, 72.20.Ht, 72.70.-|-m 
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1. Introduction 

The hopping transport of quasi-localized electrons in disordered conductors and 
semiconductors has been studied for many years - for comprehensive reviews, see 
Refs. [H El IHj. The more recent observation [H IH] that hopping transport may 
implement the quasi-continuous ("sub-electron") charge transfer, hence providing a 
possible solution to the random background charge problem in single-electronics [H], 
has renewed interest in this phenomenon, with an emphasis on the shot noise of the 
hopping current [3 |HJ IHl HHI- The objective of this paper is to present the results of 
an extension of our previous numerical studies of 2D hopping [HI HDj to the case of 
substantial Coulomb interaction of hopping electrons. Just as in the case of negligible 
interaction ITU], the use of advanced algorithms of spectral density calculation and 
modern supercomputer facilities has allowed us not only to obtain more complete and 
exact results for average characteristics of hopping transport (including the dependence 
of the dc current on temperature and electric field), but also to calculate the spectral 
density of current fluctuations at low temperatures. 

In order to explain our new findings, we have to start with a brief summary of the 
basic prior results. 

1.1. Coulomb Cap 

Most theoretical discussions of the Coulomb interaction effects on hopping are based 
on the notion of the Coulomb gap in the electron energy spectrum. Generally 
speaking, substantial Coulomb interaction makes the single-particle energy meaningless. 
However, the introduction P] of the effective single-particle energy e, which includes the 
contribution from the Coulomb interaction with other electrons, immediately leads to 
a "soft" gap in the single-particle density of states u (e) at e ~ /x, where /i is the Fermi 
level. In the case of 2D conductors with the 3D Coulomb interaction law, which is the 
focus of our current work, simple arguments j21 IHj yield 



where e is electron charge, k is the dielectric constant of the insulating environment 
and c is a dimensionless constant. Equation is valid only when the 2D density of 
states V (e) is much smaller than the "seed" density of states z/q; for larger e there is a 
continuous crossover to vq. The effective width A of the Coulomb gap can be estimated 
from the natural condition v (A) = z/q, resulting in 



CK^ 

A self-consistent approach jS] allows a more rigorous evaluation of the Coulomb gap 
width, giving c = 2/-n. 




(1) 
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1.2. DC Transport Characteristics 



At low applied electric fields -E, the average current (J) is a linear function of E, i.e. 
the 2D ("sheet") dc conductivity a{T,E,x) = (I) /EW (where W is the width of 
the conductor) is independent of E. For not very high temperatures (T <^ Tq, where 
ksTo = 1/h'oa'^ and a is the localization radius), the ratio a/ao (where uo is a constant 
characterizing the sample) depends only on two dimensionless parameters: ratio T/Tq 
and parameter x = e^vod/f^ characterizing the Coulomb interaction strength. The 
relation between these two parameters determines two possible variable-range hopping 
transport regimes. 

If the Coulomb interaction is weak (x^ <^ T/Tq), the average length r{T,E,x) of 
the so-called "critical hops", which connect percolation clusters and hence determine 
the current, may be found from the Mott theory 

1 /3 

r(T,0,0)(x a. (3) 



In this case the conductivity is 0121 E]: 



A (T, 0, 0) exp 



5(T,0,0) 



To 
T 



1/3 



(4) 



where A (T, E, x) is a dimensionless, model-dependent, slow function of its arguments, 
while B (T, E, x) is usually treated as a constant, but in general may be also a weakly 
dependent function of its arguments. 

On the other hand, if the Coulomb interaction is strong {x^ ^ T/Tq), the critical 
hops are longer [21 El : 

r (T, 0, x) oc 



xTq 



1/2 



and the dc conductivity is suppressed [2 EI : 



a 



A(T,0,x)exp 



B{T,0,x) 



xTq 
T 



1/2 



(5) 



(6) 



In the case of relatively high electric fields {E ^ Et, where Ej- = ksT/ea), the 
dc current is a highly nonlinear (exponential) function of the applied electric field E. If 
the field is not extremely high [E <^ Eq = l/euoa^), i.e. in the variable-range hopping 
domain, we can again distinguish two different transport regimes. 

If the Coulomb interaction is weak {x^ <^ E/Eo), one can neglect the effects of 
Coulomb interaction to evaluate the critical hop length 

E 

In this case, the dc conductivity is [121 El CH US HHI 



r (0,E,0) oc 



a. 



(7) 



a 



A (0, E, 0) exp 



5(0,^,0) 



E 



1/3 



(8) 
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In the opposite limit (x^ ^ E/Eq) 




(9) 



and the dc conductivity is lower pl!Hj : 



B(0.i5.x)^) 



1/2 



A (0, E, x) exp 



(10) 



1.3. Current Fluctuations 

At low temperatures ^7], the dynamical fluctuations of the current flowing through a 
mesoscopic system are more sensitive to the charge transport mechanism peculiarities 
than the average transport characteristics, and therefore may provide additional 
information about the conduction physics fHJ 120] • If we refrain from the discussion 
of the quantum fluctuations at extremely high frequencies, two basic frequency ranges 
have to be distinguished. At very low frequencies, one can expect the l//-type noise 
that is observed experimentally in a wide variety of conductors - see, e.g., Ref. |18j . 
In most cases the noise scales approximately in accordance with the phenomeno logical 
Hooge formula fHll^- For a 2D conductor, this formula can be presented as 



where Si{f) is the current spectral density, L is the length of the conductor (along 
the current flow) and C (/) is either a dimensionless constant or a weak function 
of the observation frequency /. In particular, many studies JH] have found that 
C (/) // oc 1//^, where p is typically between 1 and 2. For the particular case of hopping 
conduction, two major theories of 1// noise have been suggested, based, respectively, 
on "carrier number" fluctuations j23 and "mobility" fluctuations [23 1211 

possible origins of the noise. Unfortunately, both theories have been developed for the 
case of substantially nonvanishing temperatures, for which an accurate numerical study 
of noise is difficult even with currently available advanced simulation algorithms and 
supercomputer resources. 

At relatively high frequencies, the noise spectral density is a very slow (practically 
constant) function of /, and may be considered as a mixture of thermal fluctuations and 
shot noise. In the most interesting case of sufficiently low temperatures, the thermal 
fluctuations are negligible, and the broadband fluctuations are completely due to electric 
charge discreteness (shot noise). 

An emphasis of most recent studies has been on the suppression of the shot noise 
with respect to its Schottky value 2e (J). In particular, such suppression is a necessary 
condition for quasi-continuous charge transfer at relatively high frequencies [HE]. If 
the current spectral density Si{f) is flat at / — > 0, it may be characterized by the Fano 
factor 



Si if) a' C{f) 



(11) 



{If LW f 



F = 



5,(0) 



(12) 



2e(/) 
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so that the term "shot noise suppression" means that F < 1. Previous theoretical 
studies of shot noise at hopping in artificial (space-ordered) ID |2Zl E] and (both space- 
ordered and random) 2D |21 systems have shown that the shot noise may be, indeed, 
suppressed, obeying 



where is a scaling constant interpreted as the average percolation cluster length (i.e. 
the average distance separating critical hops |21I2H]) and a is a positive exponent. In 
fact, at T in the limit of negligible Coulomb interactions, our prior results 
show that Lc obeys the law 



where J is a dimensionless constant of the order of 1, and the value of the numerical 
exponent is /i = 0.98 ± 0.08, consistent with the estimate /x ~ 0.91 based on directed 
percolation theory (TUl 123 123 EO] . 

Considering a very long conductor, one might suspect that the electron motion in 
distant parts should not be correlated. This assumption immediately leads to a = 1 [20]. 
However, both analytical and numerical results [71I2ZI show that at ID hopping without 
Coulomb interaction, a may be as low as 1/2. This nontrivial result may be interpreted 
as a consequence of an essentially infinite correlation length in ID conductors, due to 
the on-site interaction of hopping electrons. Even more surprisingly, the exponent a 
may be substantially below 1 even in 2D conductors. For systems on a regular lattice, 
and without the Coulomb interaction, numerical modeling [H] yields a = 0.85 ± 0.02. 
In our recent work [lOj, this finding has been confirmed for 2D hopping in conductors 
with completely random distribution of localized sites both in space and in energy. Our 
most accurate result was a = 0.76 ± 0.08, i.e. significantly below 1. 

It has not been immediately clear how the inclusion of Coulomb interaction effects 
might affect this result. For ID hopping with increasing strength of the Coulomb 
interaction, numerical results [7j show a crossing over from nearly 1/2 up to 1. The 
similar behavior might be expected for 2D hopping, because the long-range correlations, 
apparently responsible for the difference between a and 1, should be suppressed by 
Coulomb interaction effects, provided that the conductor length L is larger than a 
certain crossover length determined by the interaction constant x- Unfortunately, recent 
experiments jHJ [Ml |H21 could not help in answering this question; while giving a reliable 
confirmation of the shot noise suppression in longer conductors, their accuracy is not 
sufficient to resolve a possible (relatively minor) deviation of a from 1. 

The resolution of the problem of shot noise suppression in long conductors has been 
the main motivation for the numerical experiment described in this paper. However, 
since the calculation of dc transport characteristics is computationally much less 
demanding that that of current noise, we have used this opportunity to obtain accurate 
values for the slow functions A and B for the same model of 2D hopping. 




(13) 




(14) 
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2. Model 



We have studied broad 2D rectangular conductors {W 3> L^) with "open" boundary 
conditions on the interfaces with well-conducting electrodes [HI HO]- The conductors 
are assumed to be "fully frustrated", with a large number of localized sites randomly 
distributed over the conductor area. At the sites, the corresponding electron "seed" 
eigenenergies e^^^ are also random, being uniformly distributed over a sufficiently broad 
energy band 25, so that the 2D density of states z/q is constant at all energies relevant 
for conduction. 

The carriers are permitted to hop from any site j to any other site k with the rate 
7jfc = Tj-fcexp (^-^j , (15) 

where rjk is the site separation distance and Vjk contains the energy dependence (see 
below). Such exponential dependence on the length of a hop is standard for virtually 
all theoretical studies of hopping [32] . Following our prior work P HU] , we take Eq. (fT3|l 
literally even at small distances rj^ ~ a. The energy dependence of T is given by the 
usual formula [THj 



hTjk {AUjk) = g- 



(16) 



1 - exp {-AUjk/kBT)' 

where 5^ is a dimensionless parameter which determines the 2D conductivity scale 
(To = ge'^/h while AUjk is the difference of the total system energy before and 
after the hop from site j to site k: 

AU,k = Uj -Uk + eEr.fc. (17) 

Here U is the total internal energy of the system, including the effects of Coulomb 
interaction: 

o2 



I 



2k 



^i^T + ^ ( ^« - i ) 5^ ( - ^ ) (n, ri, 



I'^l 



(18) 



where = or 1 is the occupation number of the l-th localized site. (Similar to earlier 
studies 13 EI of the Coulomb effect on hopping, we keep the system electroneutral by 
adding a background charge of e/2 to each site.) G{rj,rk) is the electrostatic Green's 
function 



Girj,rk) = 



1 



{2nL + Xk 



(19) 

'{2nL + Xk + Xjf + {yk - yjf _ 

which includes the effect of image charges representing the screening effect of external 
electrodes modeled as ideally conducting semi-spaces. 
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For practical calculations, we do not need to evaluate U from Eq. (fTH|). because this 
equation may be used to rewrite Eq. ()17p in the explicit form 

AU,, = ef - ef + eEr,, + ^ E " ^) ^ (^^' ^0 

- 7 E - ^) ^ (^^' + 7^ (^^■' • (20) 

The numerical study has been carried out by using the classical Monte Carlo 
technique based on the algorithm suggested by Bakhvalov et al. [SHI, which has become 
the de facto standard for single-electron tunneling simulations jSH]. In most cases, 
the calculated variables are averaged over several (many) conductors with independent 
random distributions of locahzed sites in space and energy, but the same macroscopic 
parameters. The spectral density of current fluctuations is calculated using the advanced 
algorithm described in detail in Ref. 

3. Results 

In order to classify the physical regimes of hopping behavior, it is useful to note that 
our model has four relevant energy scales: 

(i) 1 / v^a^ describes the energy spectrum discreteness, 

(ii) eEa is the scale of the electric field strength, 

(iii) e^/ Ka = xj v^c? characterizes the Coulomb interaction strength and 

(iv) k^T is the scale of thermal fluctuations. 

Our primary interest is in transport, especially its dependence on the applied 
electric field i?, so that instead of comparing eEa with the other three energy scales, 
we prefer to speak about three characteristic values of electric field, which should be 
compared with the actual E: 

T 1 1 1 

eaET = keT = — x -, eaEo = -, eaE^ = — = X x -. (21) 

We are not interested in the case of extremely high temperatures, so that we will always 
assume that T ^ Tq, i.e. E^ <^ Eq. On the other hand, the relative position of points 
Ec and Eq on the field axis is determined by the normalized parameter of the Coulomb 
interaction strength: 

E,/Eo = x = eVoa/K. (22) 
3.1. Coulomb Cap 

In order to understand the peculiarities of Coulomb interaction effects in our model, we 
started with a calculation of the single-particle density of states for the case of T = 
and E = 0. Indeed, in all the Coulomb gap analyses we are aware of, the electrostatic 
boundary effects have been ignored by assuming 

G{r,,rk) = —. (23) 
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On the contrary, in our Green's function (|19p the image charge contribution may be 
substantial, so it has been essential to understand how this contribution affects the 
Coulomb gap formation. 

Following the Coulomb gap literature OEl, we define the effective single-particle 
energy of an electron on site j as 



The calculations of the single-particle density of states in the ground state of a 
system (in which all /^Ujk are negative) require its "annealing" . In our case the annealing 
is facilitated by the fact that our model allows hopping between any pair of sites. This 
is why the natural relaxation of the conductor at T = and E = {] gave the results 
undistinguishable from those obtained after an explicit annealing procedure (see, e.g., 
Ref. 03)- Since the used Monte Carlo algorithm is not slowed down when all the 
transition rates are very low, the relaxation could be simulated very quickly. 

Figure [T] shows our typical results for the single-particle density of states. The soft 
Coulomb gap at sufficiently low energies is clearly visible. The effects of screening by 
the external electrodes are shown in Fig. ^a). The data labeled "Screened" correspond 
to the full Green's function (fT^ . which includes the electrostatic screening effects of 
the external electrodes, while the results for the simple approximation are marked 
"Unscreened" . The results show that for conductors of sufficiently large size, screening 
has virtually no effect on the Coulomb gap formation. In this limit, the linear part 
of the u (e) dependence corresponds to Eq. ((H) with the self-consistent equation result 
c = 2/7r ^ 0.64 cited above. 

Figure CJb) shows the single-particle density of states for three different values of an 
important technical parameter, the half-bandwidth B of the seed energy band. One can 
see that the value of B does not affect the single-particle density of states well inside the 
Coulomb gap, but may influence the results at larger energies, so that B should always 
be chosen properly in each particular case. 

All the results presented below have been obtained for conductor size L x W and 
energy bandwidth 2B so large that the effects of screening and finite number of states 
are negligible. 

3.2. DC Transport Characteristics 

Figure El shows our Monte Carlo results for the dc conductivity cr as a function of 
temperature T for two values of the Coulomb interaction strength parameter x- The 
results for X = coincide with those discussed in our previous work [10]. In particular, 
for sufficiently low temperatures {Et -C Eq), the Monte Carlo data may be fit by 
Eq. (j3)) using a simple power law in temperature T to estimate the pre-exponential 




(24) 



Note that our basic Eq. (j2(Jj) may be conveniently rewritten in terms of ej: 
AUjk = Sj - ek + eEvjk H G (r^, r^) . 



(25) 
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Figure 1. Single-particle density of states v (e) /vq averaged over a large number of 
conductors at x = 0.5, T — and E — Q for (a) several conductor lengths L for large 
width W/a — 40 at fixed half-bandwidth of the seed energy band [B ~ 1) with and 
without the screening due to electrostatic boundary effects, and for (b) several values 
of the half-bandwidth B for sufficiently large conductors with screening. The straight 
lines correspond to Eq. JQ) with a = 2/it k, 0.64. Curves are only guides for the eye. 

(model-dependent) function 74(T, 0,0) = (23.4 ± 1.3) (T/To)^°'^^^° °^'' and a constant 
for B (T, 0, 0) = 2.0 ± 0.2. (See Ref. [TUj for a detailed discussion of this result.) 

On the other hand, the results for x = 0-5 show that in the case of substantial 
Coulomb interaction, the temperature dependence of conductance at low temperatures 
(T/Tq ^ x^) follows the Efros-Shklovskii variable-range hopping result ©• The best 
fitting gives A (T, 0,0.5) = (10.7 ± 1.3) (r/To)^^-°^^°-^^^ and 5 (T, 0,0.5) = 1.4 ± 0.3. 
This is in a reasonable agreement with the following values for B (in our units): 3.1 
following from an approximate analysis based on percolation theory 4.8 found 
by evaluating an approximate integral over critical hops using a lattice model jHS] 
and 2.9 obtained for a narrower range of temperatures using numerical (Monte Carlo) 
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Figure 2. Linear dc conductivity a for negligible Coulomb interaction and finite 
Coulomb interaction as a function of temperature. Points show the Monte Carlo results 
which were obtained by the direct averaging of current calculated for a large number 
(from 20 to 96) of conductors with a random distribution of localized states, but the 
same macroscopic parameters. The sample size ranged from 20a x 14a to 80a x 50a, 
depending on x ^^id T. Thin dashed lines are only guides for the eye, while the thick 
solid lines correspond to the best fits of the data by Eqs. and 

simulations on a uniform periodic lattice with randomly distributed energies |4()j . 
(Unfortunately, the above values had no uncertainty reported.) The difference between 
our result and the reported values is probably due to the differences between details of 
the used models - see Sec. 2 above. 

For very high temperatures [Et > Eq), the exponential temperature dependence 
of variable-range hopping theory cannot give a good description of the results, because 
in this case transport is dominated by very short hops with lengths of the order of the 
localization radius. However, even at these temperatures. Coulomb interaction effects 
lead to a drop in dc conductivity. 

For higher electric fields {E > Et)-, the dc current (/) increases faster than E, 
therefore dc conductivity, defined as cr = (/) j EW , becomes nonlinear and begins to 
increase with E - see Fig. |3] (More extensive data for the case of negligible Coulomb 
interaction, x = 0? ^^e shown in Fig. 2 of Ref. pH].) 

For T — > 0, our results ^U] for the case of negligible Coulomb interaction {x = 0) 
exhibit an exponential field dependence for not very high fields {Et <^ E <^ Eq) and 
may be fit for the entire range of fields considered by Eq. (jH} using (analogous to the 
low field simple power law in field E to estimate the pre-exponential (model- 

dependent) function A{0,E,0) = (9.2 ± 0.6) (i?/i?o)*'°'^°^°'°^'' along with constant 
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Figure 3. Nonlinear dc conductivity ct as a function of electric field E for several 
values of temperature T and Coulomb interaction strength x- Points are Monte Carlo 
results averaged over a large number (from 20 to 96) of conductors of the same size 
(ranging from 20 x \Aa? to 800 x 500a^, depending on x, T and E). Solid symbols 
show results for T = 0, while open symbols correspond to T 7^ 0. Thin dashed lines 
are only guides for the eye. Thick solid lines are the fits to the T — Q results using 
Eqs. (© and 

B{0,E,0) = 0.65 ± 0.02. (See Ref. [10] for a detailed discussfon of this result.) The 
corresponding results for x = 0.1 and x = 0.5 show that increasing Coulomb interaction 
strength suppresses the nonhnear dc conductivity, just as in the low- field case. These 
results may be well fit by Eq. ^ with A (0, E, 0.1) = (2.3 ± 0.6) (E/Eo)^°-^^^°-°^^ and 
B{0,E,OA) = 0.96 ±0.05 for x = 0.1 and A{0,E,0.5) = (3.0 ± 0.4) (E/Eo)^°-^^^°-°^^ 
and B {0, E,0.5) = 1.68 ± 0.07 for x = 0.5. It is possible that any differences due to 
fitting reflect a (weak) systematic dependence on x- 

The results for very high electric fields, E > Eq (near the localization limit), do 
not obey the variable-range hopping theory, due to very short hops of the order of the 
localization radius which dominate the transport in this case. Note, however, that even 
within this range the dc conductivity decreases with increasing Coulomb interaction 
strength. 

To summarize our dc transport results, we see a very reasonable agreement 
with variable-range hopping theory within appropriate parameter ranges. Moreover, 
we believe that our supercomputer-enabled numerical modeling has given accurate 
parameters for the coefficients of these theories for our particular model. 
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Figure 4. Spectral density Si [lo) of current fluctuations at fixed Coulomb interaction 
strength x — 0-5, as a function of observation frequency uj measured in units of 
loq = g/hi'oa^, for several values of conductor length L. Each point represents 
data averaged over 48 conductor samples at fixed parameters (x = 0.5, T = and 
E/Eq = 0.07). Small points show results for W/a = 30, while open squares are for 
W/a = 60 (at L/a = 40). Thin dashed lines are only guides for the eye. 



3.3. Current Fluctuations 

3.3.1. 1// Noise Figure El shows typical results of our calculations of current noise 
at zero temperature, finite Coulomb interaction strength and fixed electric field, for 
several values of conductor length. Of particular note is that in sharp contrast with 
the negligible Coulomb interaction case we do observe a l//-type noise at / ^ 0. 
The frequency fk of the 1// noise "knee" (the crossover from this noise to a quasi- 
flat spectral density) is relatively constant (or at most grows slowly with decreasing 
conductor length). This is what could be expected from the comparison of Eqs. (fTT|) 
and (fT^ if a ~ 1. (For sufficiently large conductor width W , (/) is proportional to W , 
so that fk should also be independent of W as well, which is consistent with our results 
for different width in Fig. |31) 

In Fig. the calculation results are plotted in the form allowing their 
straightforward comparison with the Hooge scaling fHlE]- Indeed in these coordinates, 
Eq. (fTT|) with C(/)// oc 1//^ would give a straight line dropping with slope p. We 
see that our data for / are compatible with this formula, with p ~ 1.3. This 
value is close to the one calculated in Ref. jH] using a rather different model of 
hopping, apparently more adequate near the metal-insulator transition. The result 
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Figure 5. Spectral density Si (cj) of current fluctuations at T = and x = 0.5, 
normalized to the Hooge scaling factor a? (/) /LWluq, as a function of observation 
frequency uj (measured in units of luq) for several values of electric field. Each point 
represents data averaged over 48 conductor samples of the same size (ranging from 
20ax 14a to 120ax60a, depending on i?). Lines are only guides for the eye. FojiE/Eq = 
0.07, the results are plotted for a few conductor sizes, 50a x 30a, 60a x 30a and 70a x 30a 
(small points) and 40a x 60a (open squares). The results imply that the l//-type noise 
(in this normalization) is virtually size- and field-independent. 



is also compatible with recent experiments jl2j which indicate an increase of p from 
approximately 1 on the metallic side of such transition to above one on its dielectric 
side. 

Unfortunately, more accurate determination of the noise spectral density Si{f) 
for sufficiently small frequencies and/or finite temperatures has been out of our reach, 
despite the use of advanced averaging algorithms and unique supercomputer 
resources. As a result, at this stage we cannot compare our results with the existing 
theories of 1/f noise at hopping [23 123 123 123 • 

3.3.2. Fano Factor and Cluster Length If the low- frequency spectral density is constant 
(as it is for hopping at negligible Coulomb interaction [71 [HI CD]), it is naturally 
characterized by the Fano factor - see Eq. p2|) . In the presence of a l//-type noise, 
the definition of the Fano factor is less obvious. However, Figs. ID and El show that the 
fiuctuation spectrum has an exponentially broad plateau between the 1/f noise knee and 
a crossover to another, high-frequency value. (The latter crossover at higher frequencies 
exists even in the absence of the Coulomb interaction - see Ref. ^Uj for a more detailed 
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discussion.) The large length of these plateaus gives a motivation for the generalization 
of the Fano factor definition: 

(26) 

where fp is any frequency between the 1// noise knee and the high frequency crossover. 
In addition, following Ref. [T^, we may define a similar factor on the high frequency 
plateau: 

_ Si{f ^ oo) 

= 2e (/) • ^^^^ 
Figure |3] shows that neither of these factors depend on the sample width, at least for 
reasonably large W . Figure El shows the dependence of these factors on conductor length 
L. The results for F in the case of substantial Coulomb interactions agree well j33] with 
Eq. (fT!?|l with a ~ 1, in contrast with the result a 7^ 1 for negligible interaction [Hj llOj. 
The results for the high frequency counterpart agree well with the similar expression 

m 

Foo=[^\ . » (28) 



L 

where, within the accuracy of our calculations, (3 = 1. This result is similar to that for 
negligible Coulomb interaction [TO], and may be interpreted as a result of "capacitive 
division" of the discrete increments of externally-measured charge jumps resulting from 
single-electron hops through the system |17j. Figure Efb) shows that both results can 
be collapsed onto a universal scaling curve by the introduction of certain length scales: 
Lc for F and for F^o- 

In order to compare the corresponding length scales (Lc and Lh) with the proper 
measures of hop length, we have calculated the direction-weighted average [TU] along 
the field direction 

where Xjk = x^ — Xj = —Xkj is the component of the j ^ k hop length along the apphed 
field direction, and Hj^ is the number of electrons making this hop during a certain time 
interval. For not too high fields {Ft <^ F Fq) , the results in Fig. [7| for negligible 
Coulomb interaction are in good agreement with the variable-range hopping scaling 
described by Eq. ((Zj), while for substantial Coulomb interaction they follow scaling 
similar to Eq. ©. In both cases, Lh and x^mds have a similar behavior in the entire 
range of studied fields. (See Ref. [Hlj for a more detailed discussion on this result.) On 
the other hand, L^, as determined from Eq. (fT^. has a very different scaling, especially 
for lower fields {F <^ Fq) . Namely, at negligible Coulomb interaction, the results for Lc 
follow the law JHI) with J = 0.04 ±0.01 and /i = 0.98 ±0.08 ^ (see section 1.3), while 
in the case of substantial Coulomb interactions (x = 0.5), L^ also obeys Eq. (fT^ . but 
with J = 0.16 ± 0.02 and /i = 1.26 ± 0.15. 
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Figure 6. Average Fano factor F and its high-frequency counterpart Foo (Eqs. H26(l 
and (|27|) . respectively) as functions of conductor length L normalized to: (a) the 
localization length a, and (b) the scaling lengths Lc (for F) and Lh (for Foo) (see 
Fig- El below), for two values of applied field at x = 0.5, T = and W » Lc- Straight 
lines are the best fits to the data (using Eqs. H13|l and |2H1)), while dashed curves are 
only guides for the eye. 



Following the analysis of Ref. we may use the theory of directed percolation 
[2H1 12ni IHO] to predict the following scahng: 

L.^{x)( Y ", (30) 

\\{x)-Xc\J 

where {x) is approximately equal to Xrmds, is the critical hop length along the field, 
and the critical index S\\ should be close [SD] to 1.73. Due to the exponential nature of the 
variable-range hopping scaling, \ {x) — Xc\ ~ a, and depending on the regime of hopping, 
we may use the corresponding field scaling of Eqs. (j?!) or (j^l) (see Fig. [7j) to arrive at 
Eq. dm) with either fi = l{l + 6\\) ^ 0.91 or = i (l + 6\\) ^ 1.37, respectively. These 
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Figure 7. The values of parameters Lc and L/,. giving the best fitting of shot 
noise results for Eqs. and H28|) . respectively, for sufficiently large conductors 

(L, W » Lc), as functions of electric field at T = for the cases of negligible (x = 0, 
open squares and triangles) and substantial (x = 0.5, solid squares and triangles) 
Coulomb interaction. For comparison, circles show the results for the simple direction- 
weighted average hop length along the electric field direction H29(l . Dashed curves 
are only guides for the eye, while solid lines are the best fits using the variable-range 
hopping and percolation theory predictions (see the text). 



values are in a very reasonable agreement with our simulation results, thus confirming 
the interpretation of Lc as the average length of the directed percolation cluster. 

In extremely high fields {E > Eq), the length scales Lc and Xi-mds become 
comparable to one another and both approach the localization radius a. 

4. Discussion 

To summarize, we have carried out numerical simulations of 2D hopping within a broad 
range of temperature, electric field and Coulomb interaction strength. For average (dc) 
transport characteristics, our results are in general agreement with the variable-range 
hopping theories, except for the (model-dependent) cases of "ultra-high" electric field 
and/or temperature, where the hopping length becomes of the order of the localization 
radius. 

For the spectrum of current fluctuations, our results are more significant. First, 
for the case of significant Coulomb interaction we have obtained a reliable evidence of 
l//-like fluctuations, approximately obeying the Hooge scaling (fTT|) . even at T — > 0. 
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In hindsight, this result does not seem too surprising. Due to the presence of Coulomb 
interaction, random motion of the electrons during hopping transport generates a time- 
and space-varying Coulomb field, with a quasi-white spectrum, even at T = 0. The 
effect of such a randomly changing field on localized electrons aside from from the 
hopping clusters should be qualitatively similar to that of thermal fluctuations that 
may lead to the 1// noise, for example following one of the scenarios described in 
Refs. [221 1211 123 120] • Recent experiments for hopping in quasi-3D samples [IHlli^ . 
showing a very slow change of 1// noise intensity at T — > 0, seem qualitatively 
compatible with this interpretation. 

Our second important result is that in the presence of significant Coulomb 
interaction, the quasi- white noise above the 1// noise knee is suppressed according 
to the scaling law (|T3|) with a = 1 (within the accuracy of our numerical experiment). 
This result is consistent with the simple addition of mutually-independent noise voltages 
generated by (conductor) sample sections connected in series, and hence with the 
existence of a finite correlation length. On the other hand, the results ^0] for negligible 
Coulomb interactions give a = 0.76 ± 0.08 < 1, and are inconsistent with the existence 
of such length, at least on the scale of Lc. However, in both cases the constant Lc, 
participating in the scaling law (fTSj) . may be interpreted as the length between the 
critical hops, i.e. the directed percolation cluster length. 
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